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Abstract 



A convergent perturbation method using modified Lang Firsov transformation 
is developed for a two-site single-polaron system. The method is applicable for 
the entire range of the electron-phonon coupling strength from the antiadiabatic 
limit to the intermediate region of hopping. The single-electron energies, oscillator 
wave functions and correlation functions, calculated using this method, are in good 
agreement with the exact results. 
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1. Introduction 



The interaction of conduction electrons with lattice vibrations is described by 
the so called electron-phonon problem. The Holstein model Jl| is one of the funda- 
mental models which has been studied widely in this context. The model consists 
of a one-electron hopping term, Einstein phonons at each site and a site-diagonal 
interaction term which couples the electron density and ionic displacements at a 
given site. For weak electron-phonon (e-ph) coupling the frequency of the phonons 
and the effective mass of the the electron are renormalized, which are described 
by the Migdal approximation ||. For large e-ph coupling the electrons are self- 
trapped in the lattice deformation producing a small polaron. The motion of the 
electron is then accompanied by the lattice deformation. This results in a large 
effective mass or reduced effective hopping of the dressed electrons (polarons). The 
Lang-Firsov (LF) method based on the LF canonical || transformation works in 
this strong coupling region in the antiadiabatic limit. For weak coupling the Migdal 
approximation is satisfactory. However, no conventional analytical method exists 
at present which is beleived to describe a Holstein model for the entire range of the 
coupling strength. So, it would be useful to develop or identify an analytic method 
which could be applied to both the strong and weak coupling cases. Ranninger and 
Thibblin [|J made an exact diagonalization study of a two-site polaron problem 
and showed that the behavior of the polaron differs very much from that predicted 
by the classical LF method. Marsiglio || extended those calculations to the bulk 
limit in one dimension by studying the Holstein model with one electron up to 
16 site lattices. He concluded that for intermediate coupling strength neither the 
Migdal nor the small-polaron approximation is in quantitative agreement with the 
exact results. Kabanov and Ray || and Alexandrov et al. 0] noted that for t > ujq 
the adiabatic small-polaron approximation describes the ground state energy ac- 
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curately except for intermediate coupling strength. Ranninger and Thibblin M 
and Marsiglio || studied also the correlation functions using exact diagonalization 
technique with the finite size Holstein model and found that the results are non 
trivial and cannot be described by any conventional analytical method. Recently, 
de Mello and Ranninger || emphasized this point further. 

The objective of the present work is to search for an analytical method which 
may be applicable reasonably well for the major range of e-ph coupling strength. 
For our study we consider a two-site one-electron Holstein model for which exact 
results are available [|]]. Previously, we P, |m investigated the ground state energy 
and the nature of polarons in a two-site and a four-site Holstein model using 
the modified Lang Firsov (modified LF) transformation and two-phonon coherent 
states and found that the energy obtained within such method is very close to the 
exact result. In this work we develop a perturbation expansion within the modified 
LF transformation and show that this expansion converges for the entire range of 
the coupling strength from the antiadiabatic limit to the intermediate region of 
hopping. The energy and the correlation functions calculated within our approach 
are almost identical with the exact results. 



2. Formalism 

The Hamiltonian of a two-site one-electron Holstein model reads as 

H = Y en ia - Y *( C 1 -C 2< t + 4<r C la) + 9^0 Y n i°( h i + & 1) + U Y b M (1) 

i,cr <T i,<7 i 

where i =1 or 2, denotes the site. Ci a (c| CT ) is the annihilation (creation) operator 
for the electron with spin a at site % and rii a (=c\ a Ci a ) is the corresponding number 
operator, g denotes the e-ph coupling strength, hi and b\ are the annihilation 
and creation operators, respectively, for the phonons corresponding to interatomic 



vibrations at site i, ujq is the phonon frequency. In Hamiltonian (1) there is no 
spin dependent or spin reversal term so for the study of one-electron case the spin 
index is redundant. In the following we shall not use the spin index. 

Introducing new phonon operators a = (b\ + 62)/ v2 and d = (b\ — b 2 )j\[2 
the Hamiltonian is separated into two parts (H = Hd + H a ) : 

H d = ^2 era - t(c\c 2 + 4 c i) + ^o9+{ n i ~ n 2 )(d + d)) + uj Sd (2) 

i 

and 

H a = w (tfa - uj n 2 gl (3) 

where g + = g/y/2, a = a + ng + and a) = a) + ng + . 

H a describes a shifted oscillator which couples only with the total number of 
electrons n(= m + n 2 ), which is a constant of motion. The last term in Eq.(3) 
represents lowering of energy achieved through the lattice deformations of sites 1 
and 2 by the total number of electrons. 

Ha represents an effective e-ph system where phonons directly couple with 
the electronic degrees of freedom. Hd cannot be solved exactly by any analytical 
method. We now use the modified LF transformation where the lattice deforma- 



tions are treated as variational parameters R 11, 12]. For the present system 



H d = e R H d e- R (4) 

where R = A(ni — n 2 ) (eft — d), A is a variational parameter and linearly related to 
the displacement of the d oscillator. 

The transformed Hamiltonian is then obtained as 

Hd = u Q d^d + ^2 e p rii — t[c\c 2 exp(2A(d^ — d)) 

% 

+ c\c x exp(-2A(d t - d))] + ^ (g+ - A)(m - n 2 )(d + d f ) (5) 
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where 

e p = e - ^0(23+ - A) A (6) 

It may be mentioned that with an ordinary LF transformation, where one 
chooses a phonon basis of the oscillator with a fixed displacement (A = g + for this 
case), one can diagonalize the Hamiltonian in absence of hopping. The hopping 
term, containing off-diagonal matrix elements in the new phonon basis, may then 
be treated within the perturbation approach in the strong coupling and antiadi- 
abatic limit However, in order to develop a perturbation theory to be valid 
for the entire range of coupling strength one should consider a variational phonon 
basis such that the major part of the hamiltonian could be diagonalized for differ- 
ent values of the coupling strength. The modified LF transformation, where the 
phonon basis are formed by the oscillator with variable displacement, would serve 
this purpose. 

For the single polaron problem we choose the basis set (for Hj) 

\+,N) = ^=(c\ + cl)\0) e \N) 

(7) 

\-N) = -L(c\-ci)\0} e \N) 

where |+) and |— ) are the bonding and antibonding electronic states and \N) 
denotes the iVth excited oscillator state. 

Note that the last term in Eq.(5) has only off-diagonal matrix elements con- 
necting bonding and antibonding states with the change in phonon number by 
±1, 

(M, ±\co (g + - A)(m - n 2 ){d + J)\t, N) 
= (VN 5 m ,n-i + VN+1 S m ,n+i)ujo(9+ ~ A) (8) 
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while the hopping term H t = — t[c|c 2 exp(2A(d* — d)) + c^cx exp(— 2A(d* — d))] has 
both diagonal and off-diagonal elements in the chosen basis. The diagonal part of 
H t is given by, 



(N,±\H t \±,N) = T teT, 



i=0 



(9) 



JV! 



where t e = t exp(-2A 2 ) and N Cl ~- il(N _ i)r 

The diagonal part of the Hamiltonian H d (in the chosen basis) is considered 
as the unperturbed Hamiltonian (Hq) and the remaining part of the Hamiltonian 
Hi = Hd — Hq is treated perturbation. 

The unperturbed energy of the state |±, iV)is given by 



E® N = (N, ±\H \±, N) = Nu + e p Tt e 



' N 

E( 

,i=0 



(2A) 



2i 



(10) 



The general off-diagonal matrix elements of Hi between the two states |±, N) and 
|±, M) are calculated as (for (N — M) > 0) 



(N,±\Hi\±,M) = P(N,M) for even (N-M) 



(11) 



(N,±\Hi\t,M) = P(N,M) + VNu (g + -\)5 NM +i for odd (N-M).(12) 



where 



P(N,M) = T te(2\) N - M , 



M! 



M 



(N-M)\ £^ LV ; 



-M(M — 1)...(M — R + l)] 



(N - M + R)\RV 

In the following we present the perturbation corrections to the energies and the 
correlation functions for the ground and the first excited state. 
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3. The energies and the correlation functions 



A. Ground state : 

For the system considered, the state |+)|0) has the lowest unperturbed energy, 
-E'o = e p ~ te- The matrix element connecting this ground state and an excited 
state \e, N) is given by 

(2X) N 



(JV.eliZil^O) 



-t e 
-U 



(2A) JV 



+ u (g+ - A)5jv,i 



<5 6)+ for even N 



for odd N 



(13) 



The first order correction to the ground state wave function is obtained as, 



_ M9 + ~ A) - 2Xt e t e (2X) N 
^ } " [-0 + 2t e (l - 2A')] 1 ' 1} + N U,. m{E% - E®) ' + ' } 



t e (2\) 



N 



+ e ^: (0) ; F( o), i-^o (14) 



o 



where -E^^ is the unperturbed energy of the state |±, N) as given in Eq.(10). 
The second order correction to the ground state energy is given by 

F (2) = tl{2X-^{g + -X)Y " g(2AT n ,, 
[u* + 2f e (l - 2A*)] jfe,JV!(^ _£?(«») 1 J 

where e=+ or - for even and odd N, respectively. 

Now, one has to make a proper choice of A, hence choice for the displaced 
phonon basis, so that the perturbative expansion becomes convergent. In the 
usual modified LF method A is found out by minimizing the ground state energy 
of the system. Here we adopt that method and check whether it gives satisfactory 



results. From our previous studies P, [LOL 11] we know that A remains small as long 



as g + < 1, while for large values of g + (in the strong coupling limit) it approaches or 
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attains the full LF value of g + (see Table I). For small values of A the perturbation 
series involving linear Frohlich type (polaron-phonon) interaction term (oc (<?+ — A)) 
converges automatically, while that involving hopping, containing powers of 2A, 
would converge provided t < ujq or t ~ uq. For strong coupling, as A approaches to 
g + , the Frohlich (polaron-phonon) interaction term almost vanishes as well as t e 
becomes very small so the perturbation series converges. Thus, it is expected that 
the perturbation method following the modified LF transformation would work 
satisfactorily in both the weak and strong coupling limits. In this work we have 
shown that it works reasonably well for whole range of the coupling strength for 

t < UJq. 

Following the spirit of the modified LF method the value of A is found out as 
A = uJog + /(u)o + 2t e ) from minimization of the unperturbed ground state energy. 
It is interesting to note that for this particular choice of A, the coefficient of |— , 1) 
in Eq. (14) as well as the first term in the r. h. s. of Eq. (15) vanishes. In 
other words, the off-diagonal matrix element between the states |+,0) and |— , 1) 
becomes zero for the modified LF choice of phonon basis which leads to small 
perturbation correction to the ground state within this method. 

To check whether the perturbation series is converging properly we have calcu- 
lated and computed the third order correction to the energy for the ground state. 
The third order correction (Eq 3 ^) to the ground state energy is given by, 



*f = (E m {H 'X H J^ W (16) 

where the subscript k,m denote the states |±,iV) with the unperturbed energy 
E±^ N and the subscript refers to the ground state |+, 0). The off-diagonal matrix 
elements of Hi are calculated using Eqs.(ll) and (12). 

The ground state wave function of H d considering up to the second order cor- 
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rections in perturbation is given by, 

ig) = i+,o) + i4 1) ) + i4 2) ) 



which can be alternatively written as, 

\G) = \+,0)+ Yl a N \+,N)+ ]T b N \-N) (17) 

JV=2,4,.. N=l,3,.. 

The coefficients and b N are determined from Eq.(14) and the second order 
correction to the wave function. The normalized ground state wave function \G)n 
is 

\G) N = -^=\G) 

where N G is obtained as 

N G =(G\G} = 1+ £ a\+ J2 b 2 N (18) 

N=2,4,.. N=l,3.. 

Within the modified LF method the ground state wave function for the d oscillators 
is a displaced Gaussian 

<j>{x) = ^rexp[-(x-x ) 2 ] (19) 

7T4 

where, x — — (rii — n 2 )V2 A. Including the corrections due to the perturbation 
the ground state wave function for the d oscillator is obtained as, 

G{x) = G(x - x ) = (x - xo|0) + a N (x-x \N) 

N=2,A.. 

+ J2 b N (x-x \N) (20) 

JV=1,3... 

Note that G(x) and G(x) are the ground state oscillator wave functions for H d and 
H d , respectively. If the electron is located on site 1 then x = —a/2 A. 
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Correlation function calculation: 



The static correlation function (riiUi) and (ni-u 2 ) , where U\ and u 2 are the 
lattice deformations at site 1 and 2 respectively, produced by an electron at site 
1, indicates the strength of polaron induced lattice deformation and their spread. 
These correlation functions are determined as 



(niwi)o 
(niu 2 ) 



(21) 



<9 



A) "iVc. 



where 



A = {G\ ni (d + S)\G) = Yl b N \VN ajv-i + ^/NTl a N+1 



N=l,3,.. 



B. The First Excited State: 

The unperturbed energies of the states |+, 1) and | — , 0) are (e p +uj — t e (l— 4A 2 )) 
and (e p + t e ), respectively For 2t > u>o, the energy of the state |+, 1) is lower 
than that of |— ,0) for g + = 0, while it is higher for large values of g + when 
t e becomes negligible. The off-diagonal matrix element of H d between these two 
states is nonzero. Crossing of the unperturbed energies of these two states at an 
intermediate value of g + and nonzero off-diagonal matrix elements requires that 
one should follow the degenerate perturbation theory. So, linear combinations of 
the states |+, 1) and |— , 0) are formed to obtain two new elements of basis states 
so that Hd becomes diagonal in the sub-space spanned by these two states. The 
first excited state of Hd is described by one of the linear combinations which has 
lower energy. The unperturbed first excited state is given by 

|^ 0) )=a|-,0) + 6|+,l) (22) 
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The ratio (c) of the coefficients a and b and the unperturbed energy (a) of the first 
excited state may be found out from the relation 

c= a - Hll = Hl2 (23) 

where H n , H 2 2, H12 are the matrix elements of in the subspace of | — , 0) and 
|+, 1) and are given in the matrix form in the following, 





1- o> 






(o,-| 


(e P + te) 


2 At e + UJ (9+ 


-A) 




2 At e + u (g+ - A) 


uJo + e p - t e (l - 


-4A 2 ) 



(24) 



Eq.(23) gives two roots of a, the lower value of a (say aii) corresponds to the 
first excited state. 

The first order correction to the first excited state wave function is obtained 



as, 



VT+d 1 



E 



W P 



+ E 



tf=2A.. («i - £-,jv) 
W 



\-N) 



JV=3,5.. (C*l - E ( + > N ) 

where W c = t e ^(l + 2Ac - g) 
and W e = W + y/2u c(g+ - X)5 N)2 

Second order correction to the first excited state energy is given by 

\W \ 2 



(25) 



1 + c 2 



JV=2,4,-- 1^1 ~~ ^-,JvJ AT=3,5.. (,«! ~ -& +i jvJ 



+ E 



(26) 



4. Results and discussions 
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In this paper we report mainly the results of t=l.l (in a scale of cjq =1) for 
which exact results 0] are available. For the ground state the wave function and 
the energy have been calculated up to the second order and the third order in 
perturbation, respectively. For the numerical calculation we consider up to 25 
phonon states in the series of Eqs. (14), (15) and (16). It is found that except 
for very high values of g + cosi deration of 20 phonon states is more than sufficient, 
while for large values of g + (1.8-2.2) consideration of 25 phonon states is enough. 

In Table-I we have shown the unperturbed energy, the second and third order 
corrections to the ground state energy. It is seen that the magnitude of the higher 
order corrections decreases rapidly which clearly indicates the convergence of the 
series and reasonability of our approach. The second and third order perturbation 
corrections to the energy are small in both the weak and strong coupling limits and 
appreciable only in the intermediate coupling limit 1.0 < g + < 1.3 where higher 
order corrections may be necessary. 

It may be noted that for lower values of t (results for t =1.1 are shown here) the 
perturbation series converges more rapidly with smaller perturbation corrections. 
So, the present method based on modified LF transformation is expected to work 
very satisfactorily for t < u . 

In Fig. 1 we have shown the single electron energies as a function of g + for 
the ground state (calculated up to the third order) and the first excited state 
(calculated up to the second order). The results are found to be almost identical 
with the exact results by Ranninger and Thibblin (within the resolution of Fig. 1 
of Ref. 4). It should be mentioned that in a range 1.2 < g + < 1.4 the second order 
correction to the energy for the first excited state is ~ 10 — 12% of the unperturbed 
energy and so third order correction may be necessary in this region. For other 
regions the second order correction to the energy of the first excited state is small. 

In Fig. 2 we have shown the ground state wavefunction for the d oscillator as 
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a function of position x for different values of the e-ph coupling when the electron 
is located on site 1. For weak coupling (g + < 1) the wave function shows displaced 
Gaussian like single peak where the displacement is given by the modified LF value, 
Xo = —y/2 A. However, for g + =1.3 an additional prominent shoulder appears. For 
higher values of g + this shoulder takes the form of a broad peak. These results 
are completely consistent with the results obtained by Ranninger and Thibblin by 
exact diagonalization study Q. 

In Fig. 3 we have plotted the variation of the correlation functions {n\Ui)o and 
(niU2)o with g + . Our perturbation results are found to be very close to the exact 
results of Ref. j4j. It may be mentioned that the second order correction to the 
ground state wave function becomes very important in determining the shape of 
the correlation function within our method. In our method the ground state has 
no component of |— , 1) up to the first order correction, but it appears in the second 
order correction to the wavefunction. Presence of |— , 1) in the ground state has a 
significant contribution to the correlation function. The correlation functions are 
found to be very sensitive to any small correction to the wave function unlike the 
single elctron energies. We find a slight departure of our results from the exact 
results of Ref. citeRT at intermediate coupling strength. This is due to the finite 
series (up to the second order in perturbation for the wave function) that we have 
considered. 

It may be mentioned that recently de Mello and Ranninger concluded from 
their study of a two site one polaron problem that it would be very difficult for 
any analytical method to describe the Holstein model except in the extreme adia- 
batic or nonadiabatic limit. We have shown here that even for t/uj = 1.1 which 
is in between the above two limits the exact results are fairly reproducible by our 
analytical method based on the modified LF transformation and perturbation ex- 
pansion. Comparisons of the energies and the correlation functions in the ground 
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state with exact results show that this method, used here, with the second order 
perturbation in the wave function could match the exact results to a desirable 
accuracy except in a narrow range of g + from 1 to 1.3. For this narrow range of 
g+ reasonable results are obtained but one should include higher order corrections 
to obtain results to match with the exact one. It should be mentioned that as one 
decreases the value of t and move towards the antiadiabatic limit the convergence 
becomes better and the region of g + , where the perturbation corrections are appre- 
ciable, becomes narrower. In that case the perturbation method up to the second 
order correction in the wave function would describe the system for a wider range 
of g + . 



5. Conclusion 



In the present work we develop an analytical perturbation method within the 
modified LF approach to deal with an electron-phonon system for the whole range 
of e-ph coupling strength. This method is applicable from the antiadiabatic limit 
(t < uo) to the intermediate region of hopping (t ~ uj ). Considering a two-site one 
electron system we have calculated the single electron energies, the ground state 
oscillator wave functions and the correlation functions (up to the second order 
perturbation correction to the wave function) and find that the results are in the 
good agreement with the exact results. The perturbation series converges quite 
rapidly for all values of g + . 
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Table I. Variational parameter (A), unperturbed single electron energy (mea- 
sured with respect to the bare site energy e), second order correction and third 
order correction to the energy for different values of the coupling strength (g + ) for 
the ground state of two-site one polaron problem. 



9+ 


A 


*f - e 


^0 


E (3} 


0.2 


.0628 


-1.1125 


-.00007 


.00000 


0.5 


.1619 


-1.1795 


-.00275 


.00001 


0.8 


.2771 


-1.3099 


-.01767 


.00062 


1.0 


.3763 


-1.4397 


-.04564 


.00522 


1.1 


.4421 


-1.5212 


-.07266 


.01408 


1.2 


.5363 


-1.6183 


-.12168 


.03805 


1.3 


1.0202 


-1.7489 


-.28799 


.09973 


1.4 


1.3047 


-1.9874 


-.20805 


.02331 


1.7 


1.6875 


-2.8935 


-.11787 


.00094 


1.9 


1.8969 


-3.6108 


-.09082 


.00014 


2.2 


2.1997 


-4.8401 


-.06240 


.00002 
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Figure captions : 

FIG. 1. Single electron energies (in units of u = 1) as a function of the coupling 
strength (g+). Dashed curve: ground state, solid curve: first excited state. 

FIG. 2. Ground state oscillator wave function G(x) as a function of x for 
different values of the coupling strength when the electron is located on site 1. 
Solid: g + = 0.1, short-dashed: g + = 0.7, long-dashed: g + = 1.3 and dot-dashed: 

g+ = 2.0. 

FIG. 3. Plot of the correlation functions (a): (ni^i)o and (b): (niW 2 )o versus 
g + . To compare the results of Ref. (4) we use a unit of \\f^ J for the correlation 
functions. 
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